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. Abstract 

We use the vorticity formulation to study the long-time behavior of solutions 
to the Navier-Stokes equation on M^. We assume that the initial vorticity is small 
and decays algebraically at infinity. After introducing self-similar variables, we 
compute the long-time asymptotics of the rescaled vorticity equation up to second 
] order. Each term in the asymptotics is a self-similar divergence-free vector field with 

Gaussian decay at infinity, and the coefficients in the expansion can be determined 
by solving a finite system of ordinary differential equations. As a consequence of 
our results, we are able to characterize the set of solutions for which the velocity 
I field satisfies ||u(-, t)||2,2 = o(t~^/^) as t — > -|-oo. In particular, we show that these 

• solutions lie on a smooth invariant submanifold of codimension 11 in our function 

d . space. 

S 

• ^ : 1 Introduction 

X. 

■ We consider the motion of an incompressible viscous fluid filling the whole space M^. If 
no external force is applied, the velocity u(x, t) of the fluid satisfies the Navier-Stokes 
equation 

dtu + (u ■ V)u = z/Au - - Vp , div u = , (1) 

P 

where p is the density of the fluid, z/ is the kinematic viscosity, and p{x, t) is the pressure 
field. Replacing x, t, u, p with the dimensionless quantities 

X vt Lvl L'^p 

where L is an arbitrary length scale, Eq.([l|) is transformed into 

dtu + (u • V)u = Au - Vp , div u = . (2) 
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Since the length L was arbitrary, Eq.(^ is still invariant under the scaling transformation 



for some a > 0, then 



u(x, t) ^ Au(Ax, \H) , p{x, t) ^ \^p{Xx, XH) , (3) 

for any A > 0. 

As no external force is applied, it is intuitively clear that all finite-energy solutions 
of (I) should converge, as time goes to infinity, to the rest state u = 0, p = const. As 
a matter of fact, if u{x,t) is any global weak solution in L^(M^) satisfying the energy 
inequality, it is known that ||u(-, t)||j;^2 — as t — >• cx) (see Wiegner [^). Moreover, if 

i|e*M-,0)|U.<^^|^, t>0, (4) 

||u(-,t)|U. <^^^^, t>0, (5) 

where (3 = min(a,5/4). This result shows that the solutions of (0) decay to zero at the 
same rate as those of the linear heat equation, provided this rate does not exceed 
As we shall see below, the restriction /3 < 5/4 in (|^) is due to the nonlinearity in (j^) and 
to the incompressibility condition divu = 0. 

Wiegner's result raises a very natural question: can we characterize the set of solutions 
of (0) such that t^/^||u(-, t) ||j;^2 — > as t — > oo? Put differently, given a solution u{x,t) 
satisfying (|^) with (3 = 5/4:, under which conditions can we prove the corresponding lower 
bound ||u(-, t)||L2 > This problem has been intensively studied during the last 

15 years, especially by M.E. Schonbek |T^, |19], |2^, |21], who found sufficient conditions 



for such a lower bound to hold. For technical reasons, these results were established 
assuming some additional decay of the initial data uq = u(-, 0). Typically, it was assumed 
that uo e L2(R3)3 and (1 + |a;|)uo G L1(M3)3^ that (|) holds with a = 5/4. 

Very recently, T. Miyakawa and M.E. Schonbek obtained an interesting characteriza- 
tion of the "rapidly decreasing" solutions of the Navier-Stokes equation in R^, N > 2. In 
the case = 3, their result reads: 



Theorem 1.1 [|T^ Assume that uq G L2(R^)^, divuo = 0, and (1 + |a;|)uo G 



p3\3 



Let u{x,t) be a global weak solution of with initial data u(-,0) = uq, satisfying the 
bound with f3 = 5/4. For all k,e e {1,2,3}, define 



bki = / XkUe{x, 0) da; , cm = i I 

J9? Jo Jr 



Uk{x,t)ue{x,t) dx dt . (6) 



Then 

limt^/^||u(-,t)|U2=0 (7) 
if and only if there exists c > such that 

bke = and cm = cSm , k,i e {1, 2, 3} . (8) 

The proof is a direct calculation using the integral equation satisfied by the solutions of 
(0). As the authors themselves remark, this argument does not provide much intuition as 
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to the meaning of the conditions (|^). From our point of view, the most surprising feature 
of Theorem |0| is the fact that assertion (|^) is translation invariant in time, whereas 
conditions are not. More precisely, if a solution u{x,t) satisfies (|^), so does u{x,t + T) 
for any T > 0. However, the property that (1 + |x|)u be integrable is not preserved under 
the Navier-Stokes evolution (even for strong solutions, where uniqueness holds). Thus, if 
(1 + |a;|)uo e L^(R3)3, then in general (1 + |x|)u(-,r) ^ ^^(M^)^ for T > 0, so that the 
moments bke do not make sense if u{x, 0) is replaced by u{x,T). Moreover, if the matrix 
(cfcf) is scalar, this property will be lost if u(x, t) is replaced by u{x,t + T). Summarizing, 



Theorem |1.1| is a characterization of those rapidly decreasing solutions of (g) whose initial 
data lie in the noninvariant subspace {uq G L^(M^)^ | (1 + |a;|)uo G L^(M'^)^}. Nontrivial 
examples of such solutions have been recently constructed by L. Brandolese [|l|. 

In this paper, we use the vorticity formulation to study the long-time behavior of the 
solutions of the Navier-Stokes equation (H). Setting uj = rotu, Eq.(0) is transformed into 

dtuj + (u ■ V)uj - {uj ■ V)u = Au; , div a; = . (9) 
The velocity field u can be reconstructed from uj via the Biot-Savart law: 

1 /■ (£_y)A|Md,, .eM', (10) 

where A denotes the cross product in M'^. Although and are formally equivalent, 
we believe that using the vorticity formulation to compute the long-time asymptotics 
has a crucial advantage: roughly speaking, the spatial decay of uj is preserved under the 
evolution defined by (P). For instance, if (1 + IxD^u^o G L^(]R^)'^ for some m > 0, then 
(y) has a unique local solution uj{x,t^ with initial data ujq satisfying (1 + |a;|)™'Li;(-, t) G 
L^(M'^)^ whenever it exists. Again, we point out that this property does not hold for 
the velocity field u(a;, t) if m > 5/2. This is the reason why the integrability condition 
(1 + |x|)u G L"'^(R^)'^ is not preserved under evolution. 

In the sequel, we always assume that the vorticity uj{x,t^ is small and decreases 
sufficiently fast as |x| oo. The smallness assumption is not a restriction as far as the 
long-time behavior is concerned, since all global solutions of the Navier-Stokes equation 
(in the energy space) converge to zero as t oo. Moreover, this hypothesis allows to 
deal with global strong solutions of @. On the other hand, assuming that the vorticity 
decreases rapidly as |x| ^ cxd is very reasonable from a physical point of view. This is the 
case, for instance, if the initial data are created by stirring the fluid with a (finite size) 
tool. In addition, this property is very helpful to study the long-time asymptotics, since 
in parabolic equations the spatial and temporal behaviors of the solutions are intimately 
connected. 

To actually compute the asymptotics, we express the vorticity uj{x,t) in terms of 
the self-similar variables (^, t) defined by ^ = x/ y/l+t, r = log(l+t), see (|l^) below. 
Although the transformation is time- dependent, the rescaled vorticity w(^, r) still satisfies 
an autonomous equation, as a consequence of the scaling invariance (^). Linearizing this 
equation around the origin w = 0, we find that the generator A of the time evolution 
has a countable set of real, isolated eigenvalues with finite multiplicities, and that the 
essential spectrum can be pushed arbitrarily far away into the left-half plane by choosing 
the function space (i.e., the spatial decay of the vorticity) appropriately. Thus, the long- 
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time asymptotics in a neighborhood of the origin are determined, at any prescribed order, 
by a finite system of ordinary differential equations. 

This reduction procedure, or some variant of it, has been often apphed to investigate 
the long-time behavior of solutions of nonlinear parabolic or damped hyperbolic equations 
|TD[, ig, 0, [g, |2§. In the context of the Navier-Stokes equation. 



rescaling techniques were used by A. Carpio |Q, to study the vorticity equations in 
two and three dimensions. In f^, M. Cannone and F. Planchon constructed a large family 
of self-similar solutions of the three-dimensional Navier-Stokes equation. These solutions 
correspond to fixed points of our rescaled vorticity equation, but do not belong to the 
function spaces we use, because they decay too slowly as |a;| oo. In a companion paper 
13| , we follow the procedure outlined above to study the solutions of the two-dimensional 



Navier-Stokes and vorticity equations. In addition, we exploit the fact that the spectrum 
of the generator A is discrete to construct finite-dimensional invariant manifolds that are 
approached, at a prescribed rate, by all solutions in a neighborhood of the origin. 

The rest of this paper is organized as follows. In Section H, we prove the existence 
of global solutions of the vorticity equation (|) in a neighborhood of the origin, and we 
estimate their decay rate as t ^ oo. The results we obtain are comparable to those of 
M. Wiegner in |^. Sectionals devoted to the first order asymptotics. Under appropriate 
conditions, we show that 



where fi, f2, fs are explicit divergence-free vector fields with Gaussian decay at infinity, and 
61,62,^3 are real coefficients which can be computed from the initial data. Using (p!OD, a 
similar result can be obtained for the velocity field u{x, t). In Section we give a higher 
order asymptotic expansion of u:{x,t), including terms of the form (H-t)~^/^g(x/vT+T). 
This result is used in Section || to characterize the set of solutions t) of (|) for which 
the velocity field u{x,t) satisfies (|^). It is shown that these solutions lie on a smooth 
invariant manifold of finite codimension, which is tangent at the origin to a spectral 
subspace of the generator A. Intersecting this manifold with the (noninvariant) subspace 
{uq G L^(M^)^ I (1 + |x|)uo G //-"^(M^)^}, we recover exactly conditions (|]) in Theorem 
Finally, Appendix^ describes the spectral properties of the generator A, and Appendix]^ 
collects various estimates of the velocity field u in terms of the vorticity a; in weighted 
Lebesgue spaces. 

Notation. Throughout the paper, we use boldface letters for vector-valued functions, 
such as u{x,t) and u:{x,t). However, to avoid a proliferation of boldface symbols, we use 
standard italic characters for vector variables, such as x = (xi, X2, 0:3). In both cases, | ■ | 
denotes the Euclidean norm in M'^: |u| = {u1 + + u'^Y^'^, \x\ = {x^ + X2 + x'lY^'^. For 
any p G [1, 00], we denote by \ f\p the norm of a function / in the Lebesgue space Lp(M^). 
If f G LP(]R^)^, we set |f |p = | |f | \p. Weighted norms play a very important role in this 
paper. We always denote by p : M'^ — M the weight function defined by p{x) = 1 + 
For any m > 0, we set ||/|U = Ip'^/h, and ||f|U = Ip^^fb- If / G C%[0,T], Lp{R')), we 
often write f{-,t) or simply f(t) to denote the map x ^— f{x,t). Finally, we denote by C 
a generic positive constant, which may differ from place to place, even in the same chain 
of inequalities. 
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2 The Cauchy problem for the vorticity equation 

The aim of this section is to prove the existence of global solutions of the vorticity equation 
for small initial data in weighted Lebesgue spaces. We first recall a few standard estimates 
for the velocity field u in terms of the associated vorticity uj = rot u. Further estimates 
in weighted spaces can be found in Appendix p. 



Lemma 2.1 Let u be the velocity field obtained from uj via the Biot-Savart law 

(a) Assume that 1 < p < q < oo and i = i- i. If uj e L'P{M^^y, then u G L^(]R^)^, and 
there exists C > such that 

K<c\uj\p. (11) 

(b) Assume that 1 < p < 3 < q < oo, and define a G (0, 1) by the relation | = ^ + 
IfivE Lp(R3)3 n L''(R3)3, then u G L°°(R.^Y, and there exists C >0 such that 



|u|oo < . (12) 

(c) Assume that 1 < p < oo. If E L''(]R^)^, then Vu G Lp(M^)'^ and there exists C > 
such that 

|Vu|p < C\iv[p . (13) 
In addition, divu = and, z/divcj = 0, then rotu = lo. 

Proof: Part (a) is a direct consequence of (|T^) and of the Hardy-Littlewood-Sobolev 
inequality, see for instance Stein p2|. Theorem V.l. To prove (b), assume that uj ^ 0, 



and let R = {\u!\p/\u!\q)^ , where P = jz^j^ = 3/p-i ■ Using Holder's inequality, we find 

If.. 1 , If., .,1 



|u(x)| < — \u:{x -y)\ — dy + — / \u:{x -y)\ — dy 



1-a 



RT 

Finally, Vu is obtained from uj via a singular integral kernel of Calderon-Zygmund type. 



hence (|T3D follows from Theorem II. 3 in |2^. QED 
In the sequel, for any p G [1, cxd], we denote by L'''(]R'^) the function space 



3^ - {f G LP(M3)3 I div f = 0} , (14) 

equipped with the same norm as L^(R'^)^. As is well-known, the L^-norm of the velocity 
field u(x, t) is invariant under the scaling transformation (^. For the vorticity u){x,t), 
the corresponding critical space is L^/^(]R^). The following result shows that the Cauchy 
problem for (P) is globally well-posed for small initial data in L^/^(R^). 
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Theorem 2.2 There exists eo > such that, for all initial data ujq E L^/^(M^) with 
|^o|3/2 < eo, (|j has a unique solution u e C^{[0, oo),L^/^{R^)) n C0((0, oo), L°°(M3)) 
satisfying ci;(0) = ujq. Moreover, for all p G [|, +00], there exists Cp > such that 



\^it)\p< 



Cp|<^o|3/2 

t 2P 



t>0 . 



(15) 



Finally, ifu(x,t) is the velocity field obtained from uj(x,t) via the Biot-Savart law ([Tqj, 
then u(-,t) G L^(]R'^) for all q G [3, +00] ift>0, and there exists Cq > such that 



Ht)U< 



Cg|'^o|3/2 

1 __3_ 
t2 2q 



t > . 



(16) 



Proof: The proof of Theorem |2.2| foUows exactly the argument of Kato [|15| which shows 
that the Navier-Stokes equation has global solutions for small initial data in L^(M'^). 
The same argument also shows that the Cauchy problem for (P) is locally well-posed in 
L^/^(R^), without smallness assumption on the data. More generally, one can prove that 
(^) has global solutions for small data in the Morrey space M^/^(R^), see [|14 . 



Following [T^, we now introduce the "scaling variables" 

x 



T = log(l + t) 



VT+t 

If uj{x,t) is a solution of (|) and if u(x, t) is the corresponding velocity field, we set 

1 / X 



(jj{x, t) 
u(x, t) 



l + t 
1 



w 



X 



log(l+t)) , 

:,l0g(l+t) 



Then the rescaled vorticity w(^,r) satisfies the evolution equation 

= Aw — (v ■ V)w + (w ■ V)v , divw = , 
where A is the differential operator 

A = A5 + ^^V5 + 1, ^gM^. 

The rescaled velocity v is reconstructed from w via the Biot-Savart law: 



(17) 

(18) 
(19) 

(20) 
(21) 

(22) 



As in the two-dimensional case [0, we shall solve the rescaled vorticity equation in 
weighted spaces. For any m > 0, we define 



L'im) = {feL' 



< 00} 



(23) 



where 



1/2 



|P"^/|2 . 



Here and in the sequel, we denote by p the weight function p(^) = 1 + |.^|. In analogy 
with ([T^) , we introduce the space of divergence free vector fields 



L^(m) = {f G L^(m)^ I divf = 0} , (24) 

equipped with the norm ||f||m = Ip^lflb, where |f| = (/f + /| + fiY^'^- 

In Appendix ^ we show that the operator A is the generator of a strongly continuous 
semigroup e'^'^ in L^(m), for any m > 0. Since diA = (A + ^)di for i = 1,2,3 (where 
di = d^.), it is clear that die'^^ = e^e^^di for all r > 0. Thus, using the fact that 
divv = divw = 0, we can rewrite ( pO]) in integral form as follows: 

w,{t) = e^V(O) + y i d,e^^-'^'^''-i\w,{s)vi{s) - v,{s)w,{s)) ds , (25) 



(0) + V / dje^^-'^'^^-'^\wj{s)vi{s) - Vj{s)w,{s)) ds 
.7=1 



where i = 1,2, 3. The main result of this section states that, if the initial data are small, 
(p5|) has global solutions in L^(m) which decay exponentially to zero as r ^ +oo. 

Theorem 2.3 Let < yU < 1 and m > 2/i + |. There exists tq > such that, for all 
initial data Wq G L^(m) with ||wo||m < f^o, equation iWW has a unique global solution 
w G C"^([0, oo), L^(m)) satisfying w{0) = wq. /n addition, there exists Kq > 1 swc/i that 

||w(r)|U <iroe-'^1|wo||„ , r>0. (26) 

Proof: Given wq G L^(m), we shall solve (pSj) in the Banach space 



X = {w G C°([0,+cx)),L2(m)) | ||w||x = sup ||w(r) lUe'^" < oo} . 

r>0 

We first note that r t-^ e'^^wg G X; namely, there exists Ci > 1 such that 

lle^VolU < Cie-'^^llwolU , ^ > . (27) 

Indeed, if m < |, (|27|) is nothing but ( |62|) with e = m — 2fi — ^, a = 0, and n = — 1 or 0. 
If m > |, (pTl) follows from ( |63D with a = and ri = 0. 

Next, given w G C0([0, +oo), L2(m)), we define F[w] G C0([0, +cx)), L2(m)) by 

(Faw])(r) = V r9,e(--^)(^-^)(«;,(s)t;,(s) -t;,(s)«;,(s))ds , r > . (28) 
,=1 ^0 

We shall prove that F maps X into X, and that there exists C2 > such that 

||F[w]||x < Callwll^ , ||F[w] - F[w]||x < C2||w-W|U(||w|U + ||w||x) , (29) 

for all w,w G X. As is easily verified, the bounds (0), (|29| ) imply that the map w H-i> 
e'^^wo + F[w] has a unique fixed point in the ball {w G X \ ||w||x < R} if i? < (26*2)^^ 
and ||wo||m < {2Ci)^^R. Using Gronwall's lemma, it is then straightforward to show that 
this fixed point is actually the unique solution of (|25|) in the space C°([0, +cxd), L^(m)). 
Finally, since ||w||x < C'i||wo||m + C*2||w|||^ < Ci||wo||m + |||w||x, the bound ( PB| ) holds 
with Kq = 2Ci. 
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To prove (^), we use the following estimate, which is a consequence of Proposi- 
tions|/01and|X^. Assume that / : ^ M satisfies p™/ G L^/^(R^), where p(0 = 1 + 
Then there exists C3 > such that, for j = 1, 2, 3, 

||9,e-VlU<C3^^|p"^/|3/2 , r>0, (30) 

where a(r) = 1 — e^"^ and = min(p, i). Indeed, if < r < 2, then ( pO]) follows from 
(|6^) with p = 2 and q = 3/2. If r > 2, we have 

||a^.e(-i)AeVlU < Ce-^^^-i)||eVl|,n < C'e'^lp'^/ls/s , (31) 

where the second inequality is again a consequence of (j^. The first inequality in (|3ll) 
follows from ( |6^ ) with e = m — 2p — | and n = —1 if m < 3/2, and from ( |63D with n = — 1 
if m > 3/2. This proves (|0D for all r > 0. 

Given w G X and s > 0, we apply (|30D to / = Wj{s)vi{s) — Vj{s)wi{s), where 
i,j G {1,2,3} and v = {vi,V2,V3) is the velocity field obtained from w via the Biot- 
Savart law. Using Holder's inequality and Lemma p.l|, we can bound 



\p"^WjVi\3/2 < \p'^Wj\2\Vi\e < C\\Wj\\^\Wi\2 < C\\Wj\[rn\\Wi\\m , (32) 



SO that |p™'/|3/2 < C'4||w(s)||^ for some C4 > 0. Combining this bound with ( pOD and 
using (P^D, we obtain, for all r > 0, 



|i^.([w])(r)||^ < 3C3C4 



-(;.+ i)(r-s) 



< 3C3C4IIWIII / — e-^^'ds<C2\\w\\^xe-^^ 



since z/ + 1 > p. This establishes the first inequality in (p9D, and the second one can be 
proved along the same lines. The proof of Theorem |2.3| is now complete. QED 

Since the semigroup e"^^ is not analytic in L^(m), the solution w given by Theorem p.3| 
is in general not a smooth function of r. In particular, r 1— >■ w(r) ^ C^((0, +cxd), L^(m)), 
so that w is not a classical solution of (|20D in L^(m). Nevertheless, following the common 
use, we shall often refer to w as to the (mild) solution of (^) in L^(m). Remark that the 
evolution defined by ( pH ) is regularizing in the sense that w(^, r) is a smooth function of 
^ G for any r > 0. This property is well-known, and will not be proved here. We only 
quote the following result: 

Proposition 2.4 Let < p < 1, m > 2p + 1, and let w G C°([0, 00), L^(m)) be the 
solution of given by Theorem |^. 3^ . There exists i^i > such that, for allp G [2, +00], 

|p'"w(r)|p<iri(l + r->)e-'^^||wo|U , r>0 , (33) 
where p(0 = 1 + |^| and 7p = |(i - i). 

Proof: In view of (p6D , it is clearly sufficient to prove (0) for < r < 1. This can 
be done by a standard bootstrap argument, using Proposition |A.4| , Lemma ^.1| , and the 
integral equation ( pSj) satisfied by w. We omit the details. QED 
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Corollary 2.5 LetO < /i < 1 andm > 2/i + i. Letw E C°([0, oo), L^(m)) be the solution 



of ^U) given by Theorem 2.5, and let v be the corresponding velocity field. There exists 



and 



K2 > such that, for all t > 1, 

|w(.)|, < A-.e--||w„||„, /„. all ^^^^^^ f ^ t ' (34) 

|v(r)|, < A-.e--||w„||„ for all | | <;| ^ i ' (35) 

where Pm = 3^2m ^''^'^ ^™ ~ \+2m ■ Moreover, the bounds ^^), ( l^ j hold for all r > if 
p<2,q<6. 

Proof: If p > 2, (|3^ ) is a direct consequence of (pS]). If < p < 2, then |w|p < C|p™w|2 
by Holder's inequality, and (plj) again follows from (0). Using (0) and Lemma we 
easily obtain P3| ) if g > 3/2. Finally, if m > 3/2 and max(l, qm) < g < 3/2, we can pick 
n < min(m,5/2) such that q > qn- Using Holder's inequality and Proposition we 
find |v|, < C|p"v|6 < C|p"w|2 < C||w|U, and (|5|) follows. QED 



Remark 2.6 The fact that the value q = 1 is excluded in ( l^ j ifm > ^ is not a technical 
restriction. As is shown in Corollary \B. 4 , the velocity field \{^,t) is not integrable in this 



case, unless Jj^g iiWj{^, r) d^ = for all i,j G {1, 2, 3}. 

For the vorticity lj{x, t) and the velocity u(x, t) in the original variables. Corollary ^]5] 
implies, for the same values of p and q: 

|a;(t)|p<K2r'-^+i||u;o||m , |u(t)|, < 7^2^^"^+^ ||lc;o||„ , t>l. 



We now explain our motivation for introducing the scaling variables ([T7|). As is shown 
in Appendix 0, the spectrum of A acting on L^(m) can be decomposed as cr(A) = ad{A) U 
crc(A), where 

a,{A) = {-^ I A: G N*} , a,(A) = {a G C | 5?(A) < ^ - f } • 

(See Fig. 1.) Remark that the discrete spectrum crrf(A) does not depend on m, whereas 
the continuous spectrum (Jc(A) can be shifted arbitrarily far away from the origin by 
choosing m appropriately. Therefore, if m > is sufficiently large, the long-time behavior 
of the solutions of ( pOD in a neighborhood of the origin is governed by a finite system 
of ordinary differential equations. This system is obtained by projecting ( pOD onto the 
ffnite-dimensional subspace of L^(m) spanned by the eigenfunctions of A corresponding 
to the first eigenvalues Afc = — with k = 1,2, . . . , ko. 
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Fig. 1: The spectrum of the Unear operator A in L^(m), when m = 4. 

A rigorous justification of tliis reduction, using invariant manifold theory, can be found 
in [|T^ for the two-dimensional vorticity equation. Specifically, given any z/ > 0, we prove 
the existence of a finite-dimensional locally invariant manifold, which is tangent at the 
origin to the spectral subspace corresponding to the first eigenvalues of A, and which is 
approached at a rate (9(e~^'^) or faster by any solution of ( pOD that stays in a neighborhood 
of the origin for all times. This method allows, at least in principle, to compute the long- 
time asymptotics of the solutions to arbitrarily high order by studying a finite-dimensional 
dynamical system - the restriction of the rescaled vorticity equation to the manifold. 

Invariant manifolds can be constructed in the three-dimensional case also, and we use 
them in our discussion of the set of solutions of the Navier-Stokes equations which decay 
"faster than expected" in Section |^. However, for computing the asymptotics, we show 
that one can also use a different approach. Given G N* and m > fco + |, we decompose 
any solution w of ( PP] ) in L^(m) as 

fco Hk+2) 
k=l 1=1 

where aki G M and, for any k G {1, . . . , k^}, = 1, . . . , k{k+2)} is a basis of 

the eigenspace {w G L^(m) | Aw = — -^jiw}. Using this decomposition, (^) becomes 
a system of ordinary differential equations for the coefficients ake coupled to a partial 
differential equation for the remainder R. A direct analysis of this system allows to 
compute the asymptotics up to order 0{e~'^'^), where u = min(^2±2^ _ This program 
is carried out in Section |^ for k^ = 1 (first order asymptotics) and in Section ^ for fco = 2 
(second order asymptotics). 

3 First-order asymptotics 

In this section, we consider the behavior of the solutions of (^) in L^(m) with | < < |. 
In this space, the discrete spectrum of A consists of a single isolated eigenvalue Ai = — 1, 
of multiplicity 3. A convenient basis of eigenvectors is given by {fi,f2,f3}, where fj = 
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rot(G'ej). Here and in the sequel, G is the Gaussian function 

and {ei, 62, 63} denotes the canonical basis of M^. A short calculation shows that fj = piG 
for z = 1, 2, 3, where Pj(0 = A Exphcitly, 

, P^iO = l(^ j , P^iO = l(^ ■ (36) 

The vector fields pi satisfy divp^ = and rotp^ = ej. Integrating by parts, we thus find 
/ Pi ■ fj d^= rot(pi) ■ (Gej) d^ = (e^, ■ e^) / G = 6ij 




Moreover, if A* = A — ■ V) — | is the formal adjoint of A, it is easy to verify that 
A*pj = — Pi for i = 1,2, 3. The velocity fields v^' corresponding to fj are computed 
in Appendix p. In particular, we mention that |v^*(^)| ~ |^|~^ as |^| 00, so that 
v^* G L'^(M^) for all g > 1. Using these notations, any solution w of ( PD| ) in L^(m) can be 
decomposed as 

3 

w(e,r) = J]A(r)f.(0 + R(e,r), (37) 

i=l 

where 

(3^{r)= ! P.(0-w(e,r)de, 2 = 1,2,3. (38) 

Then R(-,t) belongs to the subspace Wi of L^(m) defined in (pl|), which is also the 
spectral subspace associated with the continuous spectrum {A G C | 3fJ(A) < | — y} of the 



operator A. As in the two-dimensional case fl^, the coefficients /5j obey a linear evolution 
equation: 



Lemma 3.1 Assume that m > |, and let w G C°([0, T], L^(m)) he a solution of ^C\). 
Then the coefficients (3i defined by j^SBj ) satisfy, for all t G [0,T], 

A(r) = -A(T), 2 = 1,2,3. 
Proof: Since rot(v A w) = (w ■ V)v — (v • V)w, ( pO]) is equivalent to 

= Aw + rot(v A w) , div w = . 
Differentiating (|38|) formally with respect to r and integrating by parts, we thus find 

A= / Pi-(Aw + rot(vAw))de = -A+ / e, •(vAw)de, (39) 



where we used the fact that A*pj = — pj. Since the right-hand side of (|39|) belongs to 
C°([0,T]) and depends continuously on w, the calculations above can be justified by a 
density argument. In particular, /5j G C^([0,T]) for i = 1,2,3. Finally, using the identity 
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V A w = V A rot V = ^V|vp — (v ■ V)v and the fact that divv = 0, we see that the last 
integral in ( pQ]) vanishes, hence Pi = —Pi. QED 

In particular, it follows from Lemma ^]T] that the subspace Wi is invariant under the 
evolution defined by (pOl). The remainder R in ( |37D satisfies the equation 

9^R= AR + gi((w- V)v-(v- V)w) , divR = , 

where Qi : L^(m) — > Wi is the spectral projection (for the operator A) onto the subspace 
Wi. Explicitly, 

3 „ 

Q^w = w-V'(/ p, -wdnfj. 

The following result describes the first order asymptotics of w(r) as r ^ oo. 

Theorem 3.2 Let 1 < u < ^, m > 2u + ^, and let w G C°([0, oo), L^(m)) be the solution 
of (^W given by Theorem with = 1. Then there exists K^>1 such that 



|w(r) - ^6.e-^f,|L < Kse-'^niwollm , r > 



i=l 



where hi = f^^ Pi ■ wq d^, i = 1, 2, 3. 

Proof: If w G C°([0, oo), L^(m)) is the solution of (pOD given by Theorem |2.3| and v the 
corresponding velocity field, we define Pi and R by (pTf), (pH]). By Lemma /5j(r) = 
fejC"'^ for z = 1, 2, 3. To bound the remainder R, we use the integral equation 



R(r)=e"^Ro+ / Qi$(r - s, w(s), v(s)) ds 
Jo 



(40) 



where Rq = QiWq and $(cr, w, v) = e'^'^((w ■ V)v — (v • V)w). 
By Proposition |A.3| , there exists Ci > such that 

lle^^Roll^m < Cie-'^^llwolU , r > . (41) 

Indeed, if | < m < |, ( p]) is a consequence of ( |62D with a = 0, n = 1, and e = m — 2z/ — |. 
If m > |, (^11) follows from ( |63|) with a = and n = 1. 

To estimate the integral in (^) , we proceed as in the proof of Theorem Exchanging 
V with e*^"^, we can write $ = ($i, $2, ^3), where 

3 

(!>i{a, w, v) = ^ dje'''-^-^\wjVi - vjWi) , i = 1,2,3 , 



see (ID. By (0), (H), there exists C2 > such that ||gi$(a, w, v) ||„ < Csa-^/^Hw] 
for all cr G [0, 1]. Using (gT]) and (^), we thus find 

||R(r)||^<Ci||wo|U+ / C2(r-s)-3/^(i^o||wo|U)Ms < C||wo|U , 

12 



2 

m 



for all r G [0, 1]. (Here and in the sequel, we use the fact that ||wo||m < ''^O) see Theo- 
rem If r > 1, we write R(r) = e^^Ro + Ri(r) + R2(r), where 



Ri(r) = 
R2(r) = 
Proceeding as above, we obtain 



T-l 



,(r-l- 



^)^Qi$(l,w(s),vGs))ds, 



Qi^{t - s, w(s), v(s)) ds . 



T-l 



Ri(r)|U < r Cie-'^(--i-^)C2(iroe-^||wo|U)Ms<Ce-^^|| 
Jo 

< r C72(r-s)-=^/^(Koe^^||wo|U)Ms<Ce-2^||wol 

Jt-1 



R2 r 



Wo I 



2 

m ■ 



Thus, there exists > such that ||R(r)||m < -ft'se '^^HwoHm, for all r > 0. 



QED 



Corollary 3.3 Let 1 < < | a^i-'i rn > 2z/+ 1. Lei w G C°([0, oo), L^(m)) 5e the solution 
of ^2[\ ) given by Theorem \2.3i with jj, = 1, and let v he the corresponding velocity field. 
There exists ^^^4 > such that, for all t > 1, 



w r 



1=1 

3 

v(r)-5^6,e-V' 



i=l 



< K^e ''^llwoilm , 1 < P < 00 , 

< fr4e~''^||wo||m 5 1 < 0' < 00 



(42) 
(43) 



where v^' is given by ([g^/ Moreover, the bounds ^^), hold for all r > if p < 2 , 
q<6. 



Proof: Using the analogue of Proposition |2.4| for R(^, r) and proceeding as in the proof 
of Corollary we obtain (|4^) for 1 < p < 00 and (^) ioi 1 < q < 00. To prove ( ^3l) for 
g = 1, we pick n G M such that 5/2 < n < min(m, 7/2). If is the velocity field obtained 
from R via the Biot-Savart law (|22|) , then using Holder's inequality and Proposition p.l| 
we can bound |v^|i < C|p"v^|6 < C|p"R|2 < C||R|U, and the result follows. QED 

In terms of the original variables. Corollary ^]3| shows that, for alH > 1, 
\u:{t) - u;app{t)\p < Cr^^'^+^llwolU , 1 < P < 00 , 



n{t) - Uappit)\^ < Ct -^-"^^^Wlj 



II m 



1 < g < CX3 



where uJapp{x,t), Uapp{x,t) are the self-similar vector fields defined by 



^ appy^^ 1 



X 



^app\-^j 



Remark that Uapp{-,t) ^ L^(]R'^), unless 61 = 62 = &3 = 0. 



1=1 ^ 



+ t)3/2 Vv/r+t 



X 
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4 Second-order asymptotics 



We now turn our attention to the solutions of (|20|) in L (m) with \ < m < ^. Acting on 
this space, the operator A has exactly two isolated eigenvalues: Ai = —1 (of multiplicity 
3) and A2 = — | (of multiplicity 8). Let £2 be the subspace of L^(m) spanned by the 
eigenfunctions corresponding to A2, see Appendix 0. A convenient basis of £2 is provided 
by the vector fields gj and hj^, which we now define. 

a) For ^ = 1, 2, 3, let = rot = rot(p,G) = f ((4 - leHe^ + ^iO- Explicitly, 



g2 = T 4-e?-e| , g3 = T . (44) 






Then divgj = and Agj = — fgi- By construction, the velocity field associated to gj is 
v^' = fj. In particular, v^' has a Gaussian decay as |^| 00. We also define 

qi(0 = ;^ 66 , q2(0 = o h-el ' ^3(^ = 7^ 66 | • (45) 



2 II / 2 ^ ^ 2 



66 / V 66 / V2-6' 

Then divq^ = 0, rot Qi = Pi, and A*qj = — |qj. It follows that 

/ qi ■ gj = / rot(qi) ■ fj d^= Pi ■ fj d^ = Sij . 

b) For (ij) E S = {(11), (12), (13), (22), (23)}, we define hij = d,ij + djfi. Explicitly, we 
have ha = ~C,iii for i = 1,2 and 

^ / 66 \ r I (i-il 

hi3 = jUl-6^ 1 , h23 = ^| -66 I . (46) 




-66 



3 




Then divhjj = and Ahjj = — |h.y for all [ij) G S. The velocity fields v*^*^ corresponding 
to hjj are computed in Appendix p. In particular, we remark that |v'^*^(^)| ~ |^|~^ as 
1^1 oo, so that pw^^^ ^ L^(M^)^. We also define rn = 1^1^362, r22 = — 166^1, and 



(47) 



Then divr.y = and A*rjj = — ffij for all {ij) G S. A direct calculation shows that the 
following orthogonality relations are satisfied: 






■ d^ = 5ik5ji , / ■ gfc d^ = / Qfc ■ hij d^ = . 

Using these notations, any solution w of (|20|) in L^(m) can be decomposed as 
3 3 

w(6r) = 5^A(r)f.(0 + E^*(^)g*(^)+ E 0,(r)h,,(0 + R(6r) , (48) 

i=l i=l {ij)&S 
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where /3j(r) is given by (|38| ) and 

l^{r) = [ q.(0 ■ w(e, r) , Qj{t) = [ r,,(0 ■ w(e, r) d^ • (49) 

Then R(-,t) belongs to the subspace W2 of L^(m) defined in (|6TD, which coincides with 
the spectral subspace associated with the continuous spectrum {A G C | 3?(A) < | — y} 
of the operator A. The coefficients 7^ and Qj satisfy the following evolution equations: 

Lemma 4.1 Assume that m > |, and let w G C°([0, T], L^(m)) &e a solution of 
Then the coefficients 'ji, dj defined by satisfy, for all r G [0,T], 

3 

7i(^) = , « = 1,2,3 , 

Cn(r) = -kuir) + l[ ivsi^,Tf~v,{^,Tf)d^, ^ = 1,2, (50) 

4(^) = f V,i^,T)v,i^,T)d^ , 1<2<J<3, 



where v = (ui, ^2, ^3) is the velocity field obtained from w via the Biot-Savart law (21). 



Proof: We proceed as in the proof of Lemma ^.11 . Differentiating (^Sp and integrating by 
parts, we find 

ii = [ qi-(Aw + rot(vAw))de = -^7i+ /" Pi " (^V|v|2 - (v ■ V)v) d^ 

JR3 Z J^3 Z 

= [ v-((vV)p,)de = -^7. , 

Z Jm.3 z 

because (v ■ V)pj = |ej A v ± v. Similarly, for all (ij) G S, we find 

Qj = -^(ij + f V ■ ((v ■ V) rot r,,) d^ • 
But rot ra = ^{^3^3 - ^iBi) and rot Tij = -H^iBj + ^jSi) for i ^ j, hence 

V ■ ((v ■ V) rot r,i) = ^i^l - v^) , v ■ ((v ■ V) rot Y,j) = -ViVj {% ^ j) . 

This concludes the proof. QED 
The remainder R in (^8]) satisfies the equation 

a^R = AR + g2((w-V)v-(v- V)w) , divR = , 

where Q2 '■ L^(m) — W2 is the spectral projection (for the operator A) onto the subspace 
W2, see Appendix 0. Our next result describes the second order asymptotics of w(r) as 
r — >• cx). 
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Theorem 4.2 Let ^ < u < 2, m > 2u + ^, and let w G C°([0, oo), L^(m)) be the solution 
of (2C) given by Theorem \2. Sj with /i = 1. Then there exist constants hi, Ci, dij, and 



such that ||w(r) — app{T)\\m < K^e '^^llwoHm for all r > 0, where 



i=l 



Proof: If w G C°([0, oo), L^(m)) is the solution of (|20|) given by Theorem p.3| and v is 
the corresponding velocity field, we define Pi, 7,, Qj, and R by ( 

3 

that /5j(r) = foje""^ and 7i(T) = Qe~2^^ where 



, ^), (B|). It is clear 



bi= [ P.(0 ■ wo(0 de , c, = /" q,(0 ■ wo(0 

JR3 JR3 



de, 2 = 1,2,3. 



By Corollary |2^, there exists Ci > such that |v(r)|2 < CiC ^||wo||m for all r > 0. 



Thus, it follows easily from Lemma iA. that ICiji^) — dije a"^] < ^^||wo||^ for all 



r > 0, where 



^ Jo JRS 



i = 1,2 



(52) 



0,(0) 



e^M V,{^,T)v,{^,T)d^dT , 1<2<J<3 



To bound the remainder R, we proceed exactly as in the proof of Theorem p.2| . By 
Proposition |A.3| , there exists C3 > such that ||e'^'^Ro||m < C*3e~'^^||wo||m for all r > 0. 
Using the integral equation 

R(r)=e"^Ro+ / Q2$(r - s, w(s), v(s)) ds , 
Jo 



together with the bound ||w(r)||m < Koe ^||wo||m given by Theorem O, it is easy to 
show that ||R(r)||^ < C4e"''^||wo|U for all r > 0. QED 

Corollary 4.3 Let | < z/ < 2 and m > 2z/ + |. Let w G C°([0, 00), L^(m)) be the solution 
of i^2dl ) given by Theorem \2.3{ with fi = 1, and let v be the corresponding velocity field. 
There exists Kq > such that, for all t > 1, 



|w(r) - w„pp(r)|p < iTge ''^||wo|U , 1 < p < 00 
|v(r) - Vapp(r)|g < Kae'^^Wwollrn , 1 < g < 00 , 



(53) 
(54) 



where 



Vapp(e,r) = $^6.e-V>(0 + $^Qe-^vS'(0+ ^«.e-^v'^-(0 • 



i=l 



i=l 



Moreover, the bounds ^^), 1 ^5^) hold for all t > if p <2, q < 6. 
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In terms of the original variables, Corollary ^]3| shows that, for all t > 1, 

|u(t) - Uapp{t)\g < Ct-^-^+^||a;olU , 1 < g < oo , (55) 
where uJapp{x,t), Uapp{x,t) are the self-similar vector fields given by 

1 / X \ 

UJapp{x,t) = ^^Wapp(^-^==,l0g(l + t)j , 
1 / X \ 

napp{x,t) = -^==Vapp(^-^==,log(l + t)j . (56) 

Remark 4.4 In contrast with the two-dimensional case, the second order asymptotic ex- 
pansions ofu}(x,t) and u(x,t) contain only integer powers o/(l + t)^^/^, and not resonant 
terms of the form (1 + t)~"log(l + 1). However, following [T^jJ , one can show that such 
logarithmic terms do appear in the third order asymptotics. This is the reason why the 



case u = 2 is excluded in Theorem \4.3i . In fact, if m > 9/2, the proof of Theorem ^71 
yields the estimate 



||w(r) - Wapp(r)||m < i^5(l + r)e ^||wo||m , r>0, 

which appears to he optimal. 

Two prior papers which discuss the second order asymptotics of solutions of the Navier- 
Stokes equations are those of Carpio, 0, and Fujigaki and Miyakawa, 0. As we demon- 
strate below, the results of Corollary |4.3| extend the results of these two references. The 
extension results from the fact that by imposing decay conditions on the initial velocity 
field (see the hypotheses of Theorem 0.6 in |^ and equation (1.2) of 0) certain terms 
in the approximating velocity field Mapp are forced to be zero. Thus, certain solutions 
of (0) of finite energy {i.e. of finite norm) whose asymptotics Corollary [4.3| allows 
us to compute are excluded from consideration by the decay conditions of and p. 
This is a further reason that we feel it is more natural to impose decay conditions on the 
vorticity rather than the velocity. Note that in deriving the higher-order asymptotics in 
0] (Theorem 2.2 (ii)) increasingly stringent decay conditions are imposed on the velocity 
which results in more and more terms in the asymptotics being zero. To compare the 
results of Corollary [4.3| with those of the previous references, first note that the require- 
ment that (1 + |x|)vo G L^(]R^)^ implies that that 6j = for i = 1, 2, 3 by Corollary pi . 
Furthermore, by Corollary |B.5| , the remaining coefficients in \app satisfy 



ci = 623 , C2 = 631 , C3 



n2 



and 



dll — 2 ('^33 — Cii) , (^22 — 2 ('"33 " '^22) ; di2 — —Cyi , (^23 — —C23 , C^13 — "ClS , 

where bke and Cke are defined in (^. The expressions above for dij follow easily from (|TPp, 
(0) if one remembers that 0^(0) = by Corollary |B.5 . 
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Now consider the term A(^) = Yl^j=i ^app- Since v^^ = fj = ^{ej A ^)G, 

where G(0 = (471")"^^^ exp(— |^p/4), this sum can be rewritten as 



1 3 1 



2 

where c = (01,02,03). Examining this expression component-by-component, we see that 

MO = Ig{0{c2^3 - C36) = loiOih.^s - &126) = -ihidsG + hid^G) , 

using the anti-symmetry of the hjk- The other components are treated in hke fashion and 
we find 

3 

AM) = -Y.^kAG{0 . J = 1,2,3. 

fc=i 

To treat the term B(^) = Yli{ij)<^s^i3^^^^^0)i 't^ote that using the expressions for v'^'j 
in Appendix |B], one has 

v*^- = 29,a,(V$) -I- {diG)ej + {djG)e, , 

where — A$ = G. Then using the expressions for dij from above a straightforward com- 
putation shows that 

3 

In 0, the following notation is introduced: 

^t(x) = ^G(^) , F,,,,(x,t) = iF,,,fe(^,l) , xeR^ t>0, 

where Fij^i', 1) = <9fe<9^(<9j$) + {diG)6jk. Using (|56| ) and the expressions above for A(^) 
and B(^), it is easy to verify that the velocity field Uapp in (^) satisfies 



t-l) = - ^ bkjdkEt{x) - E CkeFijk{x, t) , j = 1, 2, 3 . 



k=l k,e=l 

Thus, written in terms of this notation. Corollary [4.3| implies: 

Corollary 4.5 // in addition to the hypotheses of Corollary \4.^ , one assumes that the 
initial condition for the velocity field satisfies puQ G L-^(]R^)^, then for j = 1,2,3 one has 
the estimate 



k=l kl=l 



< Ct"^~''+^||a;o|L , t>l 



Comparing with equation (2.4) of ^ we see that this is compatible with the results of 
Fujigaki and Miyakawa. Rewriting the asymptotics in a slightly different way one finds 
that they also agree with Theorem 0.6 of 0. 
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5 The strong-stable manifold of the origin 



In this section, we assume that m > 7/2 and we consider in more detail the dynamics 
of ( pnp in the invariant subspace Wi of L^(m) defined by (pTD- If wq G Wi satisfies 
II Wo II m < ''^0! where tq > is as in Theorem p73| , the solution w(-,r) of (|20D with initial 
data Wo can be decomposed as 



w(e,r) = j;7.(r)g.(0 



i=l 



where 7^, dj are defined in 
the velocity field, we have 



J2 O,(r)h.,(0 + R(e,r) 



(57) 



9]) and R(-,t) belongs to the subspace W2 of L^(m). As for 



v(e,r) = $^7.(T)v^'(0+ E 0,(r)v^"^(O+v^(C,r) 



Setting 6j = in Theorem [4.2| , we see that ||w(-,r 



0(e- 



as r 



-00. We now 



define the local strong-stable manifold of the origin by 



W. 



loc 



Wq G L^(m) 



||wo' 



I m < '"O ' 



lim e2^ 



||$rWo||m = 0| , 



(5S 



where $rWo = w(r) is the solution of (^) in L^(m) with initial data wq. It is clear 
from Theorem p.2| that W^"'^ C Wi. On the other hand, it follows from Lemma ^]T| that 
^ W2. Indeed, the inclusion W^j"^ C W2 would mean that the integrals in the right- 



hand side of ([50|) vanish identically for vorticities w G W2, but it is easy to verify that 
this is not the case. 

By construction, Wi = W2 © V, where V is the eight- dimensional space spanned by 
the vector fields gj for i = 1,2,3 and hij for (ij) G S. Using invariant manifold theory 
as in 
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it is rather straightforward to show that lyj"^ is a smooth submanifold of Wi 
which is tangent at the origin to the subspace W2. In other words, there exists a smooth 
function / : W2 ^ V satisfying /(O) = 0, /'(O) = 0, and such that 1^]°^ = G{f) D S(ro), 
where 



Gif) = {w + /(w) I w G W2} C Wi , fi(ro) = {w G Wi I ||w|U < ro} . 

In particular, the manifold W^"'^ is of codimension 8 in Wi, hence of codimension 11 
in L^(m). By construction, W^"'^ is locally positively invariant in the following sense: 
if Wq G W!,"'^, then w(r) = $rWo G W!,"'' as long as w(r) G B{ro). If in addition 
||wo||m < ^o/Kq, where ii'o is as in Theorem 2^, then w(r) G VFj"^ for all r > and 

3 

e2'^||w(r)|L — > as r ^ +00. 



The following result is characterization of the local strong-stable manifold. 

Proposition 5.1 Fixm > 7/2, and assume thatwo G Wi C L^(m) satisfies ||wo||m < f^o, 
where ro > is as in Theorem Let w{^, r) be the solution of ( f^j wzt/i initial data 
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Fig. 2: A schematic picture of the local strong-stable manifold Wj°'^ (shaded surface). The horizontal 
plane is the infinite-dimensional subspace W2, and the vertical axis the eight-dimensional space V. The 
intersection W!,°'^ H W2, which is also infinite-dimensional, is represented by two line segments. 

Wo, and let v(,^,r) be the corresponding velocity field. Define the functions '~^i{T)^Qj{T) 
by ^4^ ) and the coefficients cm by with 



u(a;, t) 



X 



log(l + 1) 



^i+t vyiTt 

Then the following three statements are equivalent: 

1) lim e2^||w(-, r)||m = 0, namely G Wl"" . 

T— >00 

2) limt5/^|u(-,t)|2 = 0. 

t— >oo 

3) 7i(0) =72(0) = 73(0) = 0, 

Cll(O) = |(Cll - C33), C22(0) = i(c22 - C33), Cl2(0) = C12, Cl3(0) = C13, (23(0) 



(59) 



C23- 



Proof: We apply Theorem 
i = 1,2,3, so that ||w(r)||m 



with 4 < z/ < 



- ^. Since wq G Wi, we have 6, = for 
0(e~3'^) as r — > -|-oo. We shall show that statements 1), 



2), 3) in Proposition |5.1j are all equivalent to 

4) Ci = for i = 1, 2, 3 and dij = for all (ij) G S. 
Indeed, it is clear from (|5lD that ||wapp(r)||m = Ke~2'^, where i^^ = if and only if 4) 
holds. Thus 1) <S=^ 4) by Theorem [4 .21 . Similarly, it follows from (^) that |uapp(t)|2 = 
K'{1 + t)-5/4, where K' = ii and only if 4) holds. Since \u{t) - Uapp{t)\2 < C{1 + t)!"^ 
by (|55|) , we conclude that 2) <;=^ 4). Finally, using (H), (^21) and the change of variables 
(BDI), we obtain the relations 



da 

dj/i 



Cm(0) + ^(C33 - C, 



0,(0) -c 



i = 1,2 
1 < « < J < 3 . 



We also know that Cj = 7^(0) for z = 1, 2, 3. Therefore, 3) -vv- 4). 



QED 



Using Proposition [5.1| , it is easy to prove Theorem |1.1| in the case where wq = rot(uo) 
satisfies ||wo||m < for some m > 7/2. Indeed, let w(^,r) be the solution of (pO|) 
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with initial data wq, and define 7i(T), OjlT) by P9|). By Corollary [B.5| , the assumption 
puo e Li(]R3)3 implies that wq G Wi (so that w(r) G Wi for all r > 0) and that Qj{0) = 
for all (ij) G 5*. 

a) Assume first that (|^) holds, namely wq G W^"'^. Then point 3) in Proposition shows 



that Cii = C22 = C33 and Cj-,- = if i 7^ j, hence the matrix (cfc^) is scalar. In addition, 
since 7i(0) = 0, we have 6^ = by Corollary p.5| . This proves (||). 



b) Conversely, assume that (||) holds. Then 7i(0) = by Corollary [B.5| , and since Cm(0) = 
= |(cii — C33) for i = 1,2, Cjj(O) = = Cjj for i 7^ j, it follows from Proposition |5]1| that 
Wo G This concludes the argument. QED 

As is clear from this proof, if wq lies in W^"'^ and if the corresponding velocity field 
vq satisfies pvo G L^(]R^)'^, then necessarily wq G W2. Thus, from our point of view. 
Theorem [L.l| is a characterization of the noninvariant set W^"'^ fl W2, and not of 
itself. As was already observed, W^"'^ is not contained in W2, and it is not a priori obvious 
that W!,°'^ n >V2 7^ {0}! In fact, using a nice argument due to L. Brandolese, it turns out 
that iyj°'^nW2 is infinite-dimensional. Following we say that a vector field u : ^ M'^ 
is symmetric if it satisfies the following two properties: 

A) ui{xi, X2, X3) = U2{x3, xi, X2) = U3{x2, X3, Xi) for all X = {xi,X2, X3) G M^. 

B) For all i G {1, 2, 3}, Ui{xi, X2, X3) is an odd function of Xi and an even function of Xj 
for all j ^ i- 

If u is symmetric, then Au and (u ■ V)u are also symmetric. This implies, roughly 
speaking, that the space of symmetric velocity fields is invariant under the Navier-Stokes 
evolution (whenever defined). 

Assume now that v : is symmetric and that the vorticity w = rot v belongs 

to L^(m) for some m > 7/2. Then w satisfies: 
A') wi(6,6,6) = ^2(^3,^1,6) = ^^^3(6,6,^1) for all ^ G Ml 

B') For all i G {1,2, 3}, Wj(^i, ^2, ^3) is an even function of and an odd function of 
for all j ^ i- 

Using these properties together with (p^), (|49|), we discover that Pi = ji = Qj = for all 
hence w G W2. On the other hand, it is clear that the integrals in the right-hand side 
of ( pOD vanish identically if v is symmetric, so that Cki = in (^. Thus, if ||w||m < tq, it 
follows from Proposition |5.1| that w G Wl"'^. Summarizing, we have shown: 



(m) ||w||m < ''"0 ) w = rotv with v symmetric > . 

We conclude this section with a somewhat surprising observation. Let \E'j be the local 
semiflow defined by the vorticity equation (^ in L^(m) for m > 7/2. Ifa;G L^(m), let 
U^; be the velocity field obtained from lj via the Biot-Savart law (|TUp. Then 



W. 



loc 



This characterization follows from the equivalence 1) -v^ 2) in Proposition |5.1| and from 
the fact that the change of variables (|T8D , ([T9| ) reduces to the identity when t = 0. As a 
consequence, W^"'^ is locally invariant under both semifiows and although the orbits 
of the same initial point under \E'j and $^ are of course different! This curious property 
originates in the fact that in both the original and rescaled variables the strong-stable 
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manifold can be characterized in terms of the decay rate (as time goes to infinity) of 
the solutions lying in it, see for a more detailed discussion. In concrete terms, the 
observation above implies that the picture of W^"'^ in Fig. 2 is not affected at all when we 
return to the original variables using (p!8[) . 



A Spectrum of the operator A 

In Appendix A of [0, we study in detail the linear operator 



N > 1 



(60) 



acting on the function L'^{m) = {f E L^(]R^) 



< oo}, where 

1/2 



In particular, we determine exactly the spectrum of C: 



Theorem A.l |]T^ Fix m > 0, and let C be the linear operator in L'^{m), defined 
on its maximal domain. Then the spectrum of C is 



a{C) = [\ec\^{\)<^-'^]u[ 2 



ken 



Moreover, if m > ^ (^'^d if k E N satisfies A; + y < m, then fi^ 
eigenvalue of C, with multiplicity {'^^^~^) ■ 



-| is an isolated 



The eigenfunctions corresponding to the isolated eigenvalues /ifc = — | can be com- 



puted explicitly. Moreover, it is shown in |13| that C generates a Co semigroup e"^^ in 
L'^{m), and sharp estimates are obtained for the norm of e^"^ in various spectral subspaces 
of L'^(m). 

In this section, we adapt the results in [|1^ to the particular case where = 3 and 

Remark 



£ = A + 2 acts on the space of divergence free vector fields L (m) defined in 
that div(Af) = £div(f), so that A preserves the divergence free condition. The analogue 
of Theorem |A.1| is: 



Theorem A. 2 Fix m > 0, and let A be the linear operator (2J_) in L^(m), defined on its 
maximal domain. Then the spectrum of A is 



a{A) = { 



A e C 



3fJ(A) < 



m 



4 



U 



keN* 



where N* = N\{0}. Moreover, ifm > | and ifk G N* satisfies k+ '^ < m, then = 
is an isolated eigenvalue of C, with multiplicity k{k + 2). 
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Proof: We first discuss the discrete spectrum of A. Fix A; G N, and take a = (ai, ^2, as) G 
N'^ such that |a| = ai + 02 + = k. Then the Hermite function 0a : — * M defined by 

0, = = dt'drOrC , where G(0 = e"!?!'/^ , 

is an eigenfunction of C with eigenvalue — |. Let Ek = span{0Q | a G , |a| = fc} and 

£k = {i = (/i,/2,/3) I divf = , /, G Efc for z = 1,2,3} . 

By construction, C L^(m) for all m > 0, and any i E £k satisfies Af = — ^i^f. 
Moreover, using the characterization of Sk in Fourier variables (see ^^), it is not difficult 
to show that dim(£^fc) = k{k + 2). In particular, for any fc G N*, Afc = — is an eigenvalue 
of A with multiplicity (at least) k{k + 2). 

Next, fix A G C such that 3ft(A) < 1/4 and -(A + 1) ^ N. Proceeding as in |l|], it is 
easy to verify that the function : 1^'^ ^ defined in Fourier variables by 



m 

4 2 • 



satisfies Kip\ = Xtpx and d\vip\ = 0. Moreover, ipx G L^(m) if and only if 3?(A) < ^ 
This shows that a{A) D {A G C 1 3?(A) <j-f}. 

Now, fix 72 G Z and assume that m>0, m>n + |. Let 

W^ = \f e L\m) [ CfiO = for all a G with |a| < n] . 

^ Jr3 J 

In particular, Wn = L'^{m) if n < 0. We define closed subspaces V„, Wn of L^(m) by 
Vn = ®l=i£k and 

>V„ = {f G L2(m) I /, G W„ for i = 1,2,3} . (61) 

By definition, V„ = {0} and >V„ = L^(m) if n < 0. (We recall that any f G L^(m) with 
m > I satisfies Jj^g f (^) d^ = as a consequence of the divergence free condition; hence 
Wo = L^(m).) Using again the characterization of L^(m) in Fourier variables, it is easy 
to verify that L^(m) = Vn® W„. Let P„ : L^(m) — > L^(m) be the (unique) continuous 
projection satisfying range(P„) = V„, ker(P„) = >V„, and let (5„ = 1 — P„. In particular, 
P„ = and Qn = 1 for all n G Z, ri < 0. The following estimates on the semigroup 
e"^^ = eT'^l'^^^^ are proved in |jl3| (Proposition A. 2): 

Proposition A. 3 Let a^r) = 1 — e^"^, r > 0. 

(a) Fix m > 0, and take n G Z such that ^ + | <m<n + |. For all a G and all 
e > 0, there exists C > such that 

ll^'^e-VlU<^^:^e5(^-+^)||/l^^ r>0, (62) 
for all f eWn C L'^{m). 

(b) Fix n G N U { — 1}, and take m G R such that m > n + ^. For all a G and all 
e > 0, there exists C > such that 

i|5"e^VlU<^^^e-^1|/||„, r>0, (63) 
for all f eWn C L'^{m). 
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If m and n are as in part (a) of Proposition |A.3| , it follows from (^) that 



m , r>0, 



for all f G L^(m). By the Hille-Yosida theorem, this implies that (j{AQn) C {A G 
C 1 3fJ(A) < ;| — y}- On the other hand, by construction, we have a{AP„) = if n < 
and cr(AP„) = {-1; -§;...; -^} if n e W. In particular, a{APn) n a{AQn) = 0, hence 
the multiplicity of the eigenvalue A^ {k = 1,. . . ,n) is exactly k{k + 2). Finally, since 
cr(A) C a{APn) U a{AQn), we have 



a(A)c{AGc|5?(A)<i-^}u{- 



A; + 1 



keW 



This concludes the proof of Theorem [A.2 



QED 



The estimates in Proposition [A.3| can be generalized to weighted spaces with p ^ 2. 
For our purposes in this paper, the following result will be sufficient: 

Proposition A. 4 Let 1 < q < p < oo, m > and T > 0. For all a G N^, there exists 
C > such that, for all f G L'^{m), 



Ip^a^e^Vlp < 



c 



a(r) 



2 ij 2 



q ' 



< r < T 



(64) 



where p(^) = 1 + |^| . 

Proof: See [|13|, Proposition A. 5. 



B Bounds on the velocity field 

We first list a few identities that are satisfied by the vorticity w as a consequence of the 
divergence free condition. If w G L^(]R^)^, then 

/ ^.(^)d^ = for alH G {1,2,3} , (65) 

because Wi = div(^jw). This fact is not hard to prove in Fourier variables, but has been 
overlooked in some papers on the subject until recently. If in addition pw G L^(M^)^ 
(where p(^) = 1 + |^|), then 

/ {^^w,{O+^MO)d^ = for alH,jG {1,2,3}, (66) 

because ^iWj + ^jWi = div(^j^jw). Thus only three first-order moments of w (out of nine) 
are possibly nonzero. Finally, if we assume that p^w G L^(R'^)^, then 

/ (eiO^fc(O + e,efc«^.(O+efc6^.(O)de = for alH,j,fcG {1,2,3}, (67) 

because ^i^jWk + ^j^k'Wi + ^k^Wj = div(^j^j^jtw). This means that only eight second-order 
moments of w (out of eighteen) are possibly nonzero. 
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Next, we give explicit formulas for the velocity fields corresponding, via the Biot-Savart 
law ([l^), to the first eigenfunctions of the linear operator A acting on L^(m). 

1) If m > 5/2, the first eigenvalue A2 = —1 has multiplicity three. A basis of eigenfunc- 
tions is {fi,f2,f3}, where fj = PiG and Pi,P25P3 are defined in (^). The corresponding 
velocity fields are 

9,(V$) + Ge, , ^ = 1,2,3, (68) 

where 

Remark that — A$ = G, so that divv^* = and rotv^' = rot(G'ej) = fj. A direct 
calculation shows that v^' ~ 1^1^^ as 1^1 C)0, hence v*"' ^ L-'^(M^). 

2) If m > 7/2, the second eigenvalue A3 = —3/2 has multiplicity eight. A convenient basis 
of eigenfunctions is given by {gi}i=i,2,3 and {hjj}(jj)gs, where gj are defined in and 
hij in dg). (We recall that S = {(11), (12), (13), (22), (23)}.) The corresponding velocity 
fields read 

Clearly, v^^(^) has a Gaussian decay as |^| 00, whereas v'^*^(^) ~ |^|~^- In particular, 
v^*^ G L\R^), but pv'^^^- ^ L\R^). 

Now, we assume that w G L^(m) for some m > 7/2. Then w can be decomposed as 

3 3 

w(0 = 5Z Af^(0 + $^7.&(0 + ^^^•'^^^•(^) + ^(^) ' (69) 
i=l i=i fo')e5 

where the coefficients jSi are defined in (^) and 7^, Qj in (^9]). The velocity field v 
associated to w has a similar decomposition: 

3 3 

v(o = E^^^''(^) + E^^^''(^)+ E Ci,v'^"(o+v(o , (70) 

i=l i=l fo')e5 

where v is obtained from w via the Biot-Savart law (|TD|). In view of ( PB]) and (§B|), it is 
clear that /3i = for z = 1, 2, 3 if and only if 

/ ^^Wj{Od^ = for alH,jG {1,2,3}. (71) 

On the other hand, it follows from (|49|) and ( |67D that Qj = for all {ij) G 5 if and only if 



^22 = ^13 = -2Mf2 = -2Mf3 , 1 _ 2 _ A#3 _ n 



(72) 



where 



M;k= [ UkW^{Od^ = Mlj , t,j,k e {1,2,3} . 

If in addition 7j = for i = 1, 2, 3, then Mj^. = for all k G {1, 2, 3}. 

The main result of this section is the following estimate for the velocity field in terms 
of the vorticity: 
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Proposition B.l Let w G L^(m) for some m > 0, and let v be the velocity field obtained 
from w via the Biot-Savart law (ff^j. Define the coefficients (3i by ([5^ and Qj by i^4^). 
Assume that either 

1) 0<m< 3/2, or 

2) 3/2 < m < 5/2, or 

3) b/2<m< 7/2 and = for i = 1, 2, 3, or 

4) 7/2 <m< 9/2, /3i = fori = 1, 2, 3, anc? Qj = /or (ij) G 5. 
r/ien t/iere exists C > such that 

|p"^v|6 < C|p™w|2 , (73) 

where p(^) = 1 + |,^| . 



Remark. More generally, Proposition [B.l| holds for any (not necessarily divergence free) 
vector field w satisfying p^w G L^(M'^)'^, provided that either 

1) < m < 3/2, or 

2) 3/2 < m < 5/2 and (|6|) holds, or 

3) 5/2 < m < 7/2 and (|6|), (|7lD hold, or 

4) 7/2 < m < 9/2 and (H), (|71]), (|7|) hold. 

Roughly speaking, the result means that v(^) decays like |^|~™~2 as |^| ^ cx). 

The proof of Proposition [B.l| is naturally divided into four steps. In the case 1), the 
bound (|7^) is a direct consequence of (p2|) and of the following weighted Hardy-Littlewood- 
Sobolev inequality: 

Lemma B.2 //-1/2 < m < 3/2 and 

then |p™m|6 < C|p™cu|2. 

Proof: We use the dyadic decomposition 



oo 
j=0 



where Sq = G | |^| < 1} and Bj = 2^-^ < \^\ < 2^} for j G N*. Let 

Ui = ulsi and uJi = luIb,, i G N. Clearly Ui = Y^jen^ij^ where 

If K ~il — 1' follows from (|lT]) that |Ajj|e < C|co'j|2. If \i — j\ > 2, Young's inequality 
implies that jAjjle < M'l^^M2^'^\ujj\2, where 



Ml = sup ( / — ^ dr^)'^' , M2 = sup ( / — ^ de 



2/3 
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Hi > j + 2, then |^ - r/| > |^| - \r]\ > 2'-^ - 2^ > 2'~^ for all ^ E Bi, r] e Bj. 
Thus Ml < C2-2^^(5j)2/3 < C2-2(»-i) and Ms < C2-'^' ii{B,f'^ < C for some C > 
independent of i, j, hence |Ajj|6 < C2~i'^*~-'^|co'j|2. If j > i + 2, then |^ — "^1 > 2-'"^ 
for all ^ G -Bj, r/ G -Bj, and a similar calculation shows that |A.y|6 < (72^2 0'^*) 1^7^12. 
Summarizing, we have shown that 

\ui\Q <C^^Kij\ujj\2 , zGN, 

jGN 

where Kij = 2^l*^-'l^5(*^-'). Now, by definition of the sets Bi, we have \p^Ui\ < C2™|'Uj| 
and Ip'^ujjl > C2"'^ujj\ for all ij G N. It follows that 

jeN 

where K^J^^ = 2~l*"-'l+('""3)(*--''). In particular, l-K'i™^! < 2~°'*~-^l for some a > 0, hence 
^{m) (iggnes a bounded linear operator from £^(N) into £^(N). This concludes the proof. 

QED 

To prove ( [73| ) in the remaining cases, we also need the following variant of Lemma [B.2| : 
Lemma B.3 // — ^ < m < ^ and 



Proof: We use the same notations as in the preceding proof. If |i — j| < 1, it fol- 
lows from the Hardy-Littlewood-Sobolev inequality that \Aij\ < C\ujj\Q/5. By Holder, 
l^^ile/s ^ C p{BjY^^\ujj\2 < C2^ujj\2, hence |Ajj|6 < C2^ujj\2. If \i — j\ > 2, then 
Wjle < N'l^'^N2^'^\ujj\2, where 

/ f 1 \2/3 / /• 1 n2/3 

Proceeding as above, we deduce that \Aij\Q < C2^^\''^^\2^ijj\2 for all i,j G N. It follows 
that 

where -K'j-™'' = 2~^l*~-'l"'"'^'^*~-'^. Thus, if \m\ < 1/2, K^^^ defines a bounded linear operator 
from into ^^(N). QED 



We are now ready to prove Proposition p.l| in the case 2). If3/2<m<5/2 and 
w G L^(m), then w is integrable and (^) holds. As a consequence, we can rewrite (|22|) 
in the form 
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where eijk = sign(cr) if (ijk) is a permutation a of (123), and eijk = otherwise. Using 
the identity 

we obtain 

< m-v\\e\v\ + i^-vfmi) ■ 

Thus, it follows from (|7|) that |v(0| < + c/(0)> where 

\v\\^iv)\ . 1 f \V\\^{V)\ 



drj 



Since we already know that |v|6 < C|w|2, it is sufficient to bound |p™v|6, where po(0 = 
|^|l{|^l>i}. Applying Lemma |B.2| with u{^) = cind uj^t]) = \ri\\w{ri)\, we find 

IPcT/le < \p™'^^u\q < C\p"^~^uj\2 < C|p'"w|2. Similarly, applying Lemma p.3| with m(^) = 
kl'^giO and cj(r/) = |r/||w(r/)|, we find |p™^|6 < C\p"'-'^u\6 < C\p'^-^uj\2 < C|p™w|2. 
Summarizing, we have shown that |p™v|6 < C|p'"w|2 if 3/2 < m < 5/2. 

We next assume that 5/2 < m < 7/2 and that /3j = for i = 1, 2, 3. As was already 
observed, this implies that ( [7TD holds. Thus (^) can be written in the form 

To bound the right-hand side, we observe that 

+m-vmm-v,)-\^-v\%) ■ 

We claim that 

lleP-|e-r/|^-3|e|(e-r/)|<C|r/P(|e| + |r/|). (78) 
Indeed, this follows from the identity 

\^\' - \^-v\' - mi^ ■ v) = - ■ v) - \^-v\\v\' 

- \mm-v\ + {^-v)-m , 



and from the bound 



o<\m-v\ + i^-v)-\e<l\v\\ (79) 



which is easily proved by setting = — 

Using (|77D together with (^), (|78|), we obtain 

<c\^-v\mvm\ + \v\) 
< m-v\mv\' + \^-v\'mv\') 
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Thus, it follows from (|7|) that |v(0| < + where 

^^^^ lePV le-r^P ^' ^^^^ |ePyM3 le-r/l 

Applying Lemmas p3.2| and |B.3| as in the previous case, we easily obtain \p^f\e + IPo'dle < 
C|p'"w|2. This prove ( ffSD in the case 3). 

Finally, we assume that 7/2 < m < 9/2 and that /9i = for z = 1,2,3, Qj = 
for (ij) G S*. As we already remarked, this implies that ( [72t ) holds. Using (^2]), it is 
straightforward (but somewhat tedious) to verify that 

3 „ 

J2 / {m\^ ■ v)Vj + - 5(e • vy))wk{v) dv = 0, zeN. (80) 

Assuming (pO|), we can rewrite (^) in the form 

where 

^.(e,^) = ier(o-^.) - iene-r]i'(o-^.) 

We claim that 

< c\^-m'\vm + \v\)< m-ri\m^\' + \^-^\W'\v\') ■ (82) 

Assuming ( p2D for the moment, we deduce from ( |8TD that |v(^)| < C(/(^) + g{0), where 

ff.^ = l_ [ HVMd. 9(0 = — [ MVMdr? 
^^^^ \^\'L \^-v\' ^' ^^^^ lel^Aa le-r/l 

Applying Lemmas [B.2| and |B.3| again, we obtain |p™/|6 + IPcTfi'le ^ C|p™w|2. This proves 
(^) in the case 4). 

It remains to establish (P^. We first remark that Aj = {^j—r]j)\^\'^B + Cj, where 



B = \^\'-m^-vf-m-v\%^-v)-l\^\m-v?-mv\') , 
Cj = l{m^j-v,)-^M-v\')m-vr-mv\') ■ 

In view of (^), it is clear that \Cj\ < C\^-7]\\^f\r]f{\^\ + \r]\). Thus, it is sufficient to 
show that 

\B\<Cmvm\ + \v\) ■ (83) 

To this end, we remark that 
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+ {\^-v\ - lel)(4(e ■ V? - mv\') - '^m^ ■ VM-V\ , (84) 

where 

= liim-vi + {^-v)- \^\'mm-v\ - lei) - (e ■ v)) ■ 

Using (frgp, we find \D\ < ||^||?7p- Inserting this bound into ( P^ ) and using ([78|), we 



obtain (pSf). This concludes the proof of Proposition |B.1| . QED 
We conclude with two corollaries which are used in the preceding sections. 

Corollary B.4 Assume that w G L^(m) for some m > |, and denote by v the velocity 
field obtained from w via the Biot-Savart law (f/^j. Then v G L^(R^)^ if and only if f3i = 
for i = 1,2, 3, and in this case f^^ d^ = for i = 1,2, 3. 

Proof: Without loss of generality, we assume that 5/2 < m < 7/2. For any w G L^(m), 
we have the decomposition 

3 3 
i=l 1=1 

where the coefficients are defined in (^). Then the remainder w G L^(m) fulfills the 
moment conditions (ffll). By Proposition p.l| , the corresponding velocity field satisfies 



p^v G L^(R^)^, hence v G L^(]R^). On the other hand, it is easy to verify that ^ /3jV^^ G 
L1(M3) if and only if a = Ofori = 1,2,3. Thus v G L\R^) if and only if /5i = jS^ = = 0. 
In this case, /j^g d^ = for i = 1, 2, 3 because div v = 0, see (|65|). QED 

Corollary B.5 Assume that w G L^(m) /or some m > |, and denote by v t/ie velocity 
field obtained from w via the Biot-Savart law (|7^. T/ien pv G L^(R'^)^ and only if 
Pi = for i = 1,2,3 and Qj = for (ij) G S. In this case, the matrix [hkn) defined by 

I ikV.iO^i, k,l e {1,2,3} , (85) 
is skew-symmetric, and 612 = 73, &23 = 7i; &31 = 72- 

Proof: Without loss of generality, we assume that 7/2 < m < 9/2. If w G L2(m), then 
w and V can be decomposed according to (|69|) , (f70|). By construction, the remainder w G 
L^(m) has vanishing first-order and second-order moments. Applying Proposition |B.1| , we 
deduce that the corresponding velocity field satisfies p'^v G L®(]R'^)'^, hence pv G L^(]R'^)^. 
Therefore, using the expressions above of fj, gj, and hjj, it is not difficult to show that 
pv G L^(R^)^ if and only if /?j = for z = 1, 2, 3 and Qj = for (ij) G S. In this case, the 
matrix (bke) defined by ( ^5] ) is skew-symmetric because divv = 0, see (|6^). 

Assume now that pv G L^(R^)^, and consider the vector field A(^) = q3(0 A v(^), 
where qs is defined in (^3]). Then A G L'^/^(R'^)^ and div A = ps ■ v — qs ■ w, where 
Ps = rot qs is defined in ( p6D . It follows that 

= / div A d^ = / (ps ■ V - qg ■ w) d^ = 612 - 73 , 

in view of (PB]), (^), and (^). The relations 623 = 7i and 631 = 72 are proved in a similar 
way. QED 
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